Transport coefficients from the Boson Uehling-Uhlenbeck Equation 



Erich D. GustE and L. E. Reichf] 

The Center for Complex Quantum Systems, The University of Texas at Austin, Austin, Texas 78712 

(Dated: August 7, 2012) 

We derive microscopic expressions for the bulk viscosity, shear viscosity and thermal conductivity 
of a quantum degenerate Bose gas above To, the critical temperature for Bose-Einstein condensation. 
The gas interacts via a contact potential and is described by the Uehling-Uhlenbeck equation. 
To derive the transport coefficients, we use Rayleigh-Schrodinger perturbation theory rather than 
the Chapman-Enskog approach. This approach illuminates the link between transport coefficients 
and eigenvalues of the collision operator. We find that a method of summing the second order 
contributions using the fact that the relaxation rates have a known limit improves the accuracy of 
the computations. We numerically compute the shear viscosity and thermal conductivity for any 
boson gas that interacts via a contact potential. We find that the bulk viscosity remains identically 
zero as it is for the classical case. 

PACS numbers: 67.10. Jn,51.10.+y,51.20.+d 



I. INTRODUCTION 

Kinetic equations provide a means to derive micro- 
scopic expressions for the transport coefficients appearing 
in the equations of fluid hydrodynamics. The transport 
coefficients for dilute gases at high temperature can be 
computed using the Boltzmann equation. However, when 
the temperature is lowered enough that quantum degen- 
eracy begins to affect the behavior of the gas, one must 
use the Uehling-Uhlenbeck equation jl|, Q ■ The Uehling- 
Uhlenbeck (U-U) equation is a semiclassical extension of 
the Boltzmann equation that accounts for the identity of 
the particles. 

There are historically two main approaches to com- 
puting transport coefficients. Most authors 0-3 use the 
methods outlined by Chapman and Enskog 6j. How- 
ever, there is another method 0] which allows more di- 
rect access to the microscopic hydrodynamic modes of 
the system and produces a more direct relation between 
the transport coefficients and the collision operator of the 
linearized U-U equation. We use this second method to 
compute the transport coefficients of dilute degenerate 
boson gases. 

The transport coefficients are macroscopic quantities 
which can be measured experimentally. A connection be- 
tween the transport coefficients and the microscopic in- 
teractions between particles allows one to measure atomic 
properties with a macroscopic apparatus. Conversely, 
one can make predictions on the macroscopic behavior 
of the gas based on a microscopic model. These types of 
relations are important in atomic physics, where control- 
ling the behavior of trapped cold atoms is essential. 

Transport coefficients determine the rate of entropy 
production in fluids. They are primarily responsible for 
the damping of hydrodynamic modes, and are intrinsi- 



cally linked to the relaxation of the gas to global equilib- 
rium. The non-zero eigenvalues of the collision operator 
are also linked to the relaxation of the gas to equilibrium 
[1, IU , and a direct relation exists between the transport 
coefficients and the eigenvalues of the collision operator. 

Recently, there has been much interest in the viscos- 
ity of various high-energy and exotic systems 0, [lol - [l3| . 
This is chiefly due to importance of modeling the after- 
math of high-energy collisions and the recent prediction 
of a universal limit on the ratio of viscosity and entropy 
density. Our method of calculation may be applicable to 
these system with the appropriate modifications. 

In Section|TT]we introduce the U-U equation and in Sec- 
tion lllll we linearize it and introduce an abstract velocity- 
space basis that will be used for calculation. Section IIVI 
gives a brief discussion of the hydrodynamic modes and 
their frequencies, and how these differ for the Bose gas. 
In Section [V] we derive microscopic expressions for the 
hydrodynamic frequencies by relating them the to the 
collision operator. In I VII we combine the results of the 
previous two seconds to derive explicit expressions for the 
bulk viscosity, thermal conductivity and shear viscosity 
of the gas. Section [VIII contains the result and discussion 
of our numerical calculation of the transport coefficients. 
We conclude by summarizing our results in Sec. IVIIII 



II. UEHLING-UHLENBECK EQUATION 

We consider a dilute gas of non-condensed bosons that 
interact via a contact potential U(r) = C/ <5 3 (r), where 
Uo = Anh 2 a/m, h is Planck's constant, a is the s-wave 
scattering length, and m is the mass of a particle. The 
derivation of a kinetic equation for this system is dis- 
cussed in Ref. [lj|. Above the critical temperature for 
Bose-Einstein condensation, the dynamics of the gas is 
governed by the Uehling-Uhlenbeck (U-U) kinetic equa- 
tion which can be written 
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where pi and r are momentum and position of the 
bosons, respectively, and fi — /(r, pi,i) is the phase 
space number density of bosons in the phase space vol- 
ume pi->pi+dpi and r— >r+dr at time t. It is normalized 
so that 

dpi 



and fi = /(r, ci,t). The overall rate constant is 
8ma 2 z(k B T) 2 
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f(r, Pl ,t) = N, 
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where N is the average number of particles in the gas. 

The quantity C[fi] is the collision integral and is de- 
fined 



x S 3 (pi 4 
x [/i/ 2 (l 



dp 2 dp 3 dp 4 

P2 - P3 - Pi)S(ei + e 2 - e 3 ~ £4) 
+ /3)(l + /4)-(l + /i)(l + / 2 )/ 3 /4] 

(3) 



where t\ — p\/ (2m). 

The collision integral explicitly conserves particle num- 
ber, momentum and energy. This can be seen from the 
five integrals / dp C[f(r, p, t)] = 0, / dppC[/(r,p,t)] = 
and J dpp 2 C[f{r,p,t)} = 0. 

The stationary state (long-time global equilibrium) so- 
lution of the U-U equation is the Bose-Einstein distribu- 
tion 
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where T is the equilibrium temperature in Kelvin, /i is the 
equilibrium chemical potential, and k B is Boltzmann's 
constant. Integrating / (pi) over momentum gives the 
equilibrium particle density uq of the gas 



which now depends on the fugacity and equilibrium tem- 
perature. 



III. LINEARIZED U-U EQUATION 

In computing transport coefficients, it is sufficient to 
consider the relaxation of the gas when it is close to equi- 
librium. In that case, the distribution function /(r, ci, t) 
will be a slowly varying function of r, and will be close 
to its equilibrium value. We can then linearize the U-U 
equation for small, spatially uniform deviations from 
equilibrium and determine the characteristic rates and 
modes of relaxation. This calculation was carried out 
in Ref. ||. We now focus on the transport coefficients 
which require that we deal with spatially varying devi- 
ations from equilibrium. We linearize the U-U equation 
by writing 

/(r,c 1 ,t)-/ 1 ° + / 1 (l + / 1 °)<i>(r,c 1 ,t), (10) 

where /f = z/(e Cl — z) and $(r, Ci,t) contains infor- 
mation about the small deviations from equilibrium and 
satisfies |$| < 1. We substitute (jTUJ) into flSJ) and neglect 
terms of quadratic and higher order in $ and obtain 
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f°(pi) = -3-Li 3/2 (2), (5) where $1 = $(r,ci,i) and 



2ttH 2 



is the thermal wavelength, 



e n/(k B T) j g f U g ac ity and Li 3 / 2 (z) is a polylogarithm. 
Polylogarithms appear repeatedly for degenerate gases 
and are defined by 
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For simplicity we will use the notation tr„ = Li »+i (z). 

It is useful to introduce a dimensionless momentum 
c i = Pi/{mvT) where v T — y / 2k B T/m. Then the 
Uehling-Uhlenbeck equation can be written 
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(12) 

We now take the Fourier transform of both the space 
and time dependence of (jTTJ) according to 
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where the dimensionless collision integral C'[/i] is given 
by 
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where A is a constant with units of volume x time that 
makes 0(k, Ci,u>) dimensionless. Since Eq. (1111) is linear, 
the value of A is irrelevant in the following analysis. 

In terms of </>(k, ci, ui), the linearized U-U equation 
then takes the form 
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where <f>\ — 0(k, ci,lj). Equation (|14[) governs the relax- 
ation of deviations from equilibrium with wave vector k 
and frequency u). 
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A. Momentum Basis 

Since (| 14[) is linear, deviations </>(k, Ci,w) will not be 
coupled to those of different wavevectors or frequencies 
and we may suppress the dependence of <p(k, Ci,u>) on 
k and w. We can greatly simplify the notation in the 
following calculations by introducing an abstract vector 
\4>) and making the interpretation 0(ci) = (ci |<^>) where 
|ci) represents an abstract momentum basis. We define 
the inner product between two abstract vectors as 



(x\<P) = J dciw(ci)x*(ci)0(ci) 

where the weighting factor w(ci) is defined 
f 
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With this weighting factor, the scalar product is normal- 
ized to one so (1|1) = 1. Note that with this definition 
of the scalar product, (ci|c2) = w ^ Cl ^ <5 3 (ci — C2). 

If we interpret C[4>] as (ci\C\<fi), where C is an ab- 
stract operator corresponding to the collision operator, 
we can show that this definition of the weighting func- 
tion makes the collision operator symmetric in the sense 
that <ci|(?|ca) = {c 2 \C\cx). 

In the momentum basis, the collision operator C can 
in general be expressed as 

C= \ dci I dc 5 iu(ci)w(c5)|ci)C(ci,C5)(c 5 | (17) 



where C(ci,C5) — (ci\C\cc,) is the "kernel function" of 
C . With this form of the collision operator, we see that 

C[(j>} = (ci\C\4>) = f dc 2 w(c 5 )C( Cl , c 5 )0(c 5 ). (18) 



By Comparing this expression with (|12j) . we can deduce 
that the kernel function of the collision operator is 
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The kernel function is symmetric under interchange of 
Ci and C5, meaning that the operator C is hermitian and 
has real eigenvalues. 



B. Angle Basis 

We will find it convenient to introduce the "angle" ba- 
sis, |c, l,m), which is defined by the inner product with 



the momentum basis as 
1 



(ci|c, Z,m) 



-.S( Cl -c)Yr(0i^i). (20) 



ci^w(ci) 

Given this definition, we can deduce that 

(ci,h,mi\c 2 ,l2, m 2 ) = 8(ci - c 2 )Si u i 2 5 mum , 2 (21) 
and that the identity operator is 
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where in the summation, I runs from zero to infinity and 
m runs from — I to I. The angle basis |c, I, m) will provide 
an elegant and economical way to perform calculations 
that would be very tedious in the momentum basis. 

The collision operator can also be expressed in the an- 
gle basis as 
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where C£'™*(ci, c 2 ) = (ci, h, mi\C\c2, h, m 2 ) is the an- 
gular kernel function. Since the kernel function C(ci, c 2 ) 
possesses rotational symmetry, the angular kernel func- 
tion is diagonal in I and m and independent of m [9] . This 
special result allows us to write the collision operator as 
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C = 2tt J dc\ J dc 2 |ci, /, ', 
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C(ci,c 2 ) = J d(ci ■ c 2 )C(ci,c 2 )Pi(ci • c 2 ). 

-1 

and Pi is a Legendre polynomial. 

C. Spectrum of the Collision Operator 



The spectrum of the collision operator plays a funda- 
mental role in both the derivation of the hydrodynamic 
equations and the calculation of transport coefficients. 
By inspection, we can see that the collision operator has 
five zero eigenvalues which correspond to the five con- 
served quantities (particle number, momentum and en- 
ergy) in a two-body collision. The non-zero eigenval- 
ues are all positive, and approach a finite limiting value 
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The representation of C in the \c,l,m) basis shows that 
the eigenfunctions of the collision operator are states of 
definite I and m. This allows us to use the indices n, I, 
m to label the eigenfunctions \4> n ,i,m) of C and write 



(c, I', m'\(pn,l,m) = (f>n,l(c')Sl t irSn 



(26) 



where <p n ,i{c) is the radial part of the eigenfunction. This 
is a property of the collision operator that follows from 
its rotational invariance. If we define <j) n j im (c,0,(p) — 
(c\(f> n ,i,m) ■, where |c) = |c, 9, (p) in spherical coordinates, 
then the radial part of the eigenfunction is given by 



<M<0 = cV^) / dnYi* m (e,ip)ii> n ,i, m (c,e,ip). (27) 



For example, the particle number conservation eigenfunc- 
tion is 0o,o,o ( c ) = 1 an d the radial part is 4>q.o(c) = 
cy/ A-kw{c). 

The radial part of the eigenfunctions satisfy the or- 
thogonality condition 
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Finally, we express the linearized U-U equation as an 
operator equation. Without any loss of generality, we 
can assume that the wave vector k = ke z , where e 2 is a 
unit vector along the z-direction. The U-U equation then 
takes the form 

-iu\ct>)+iv T kc z \<t>) =-jC\(f>), (29) 

where c z is now an operator defined by 



c, c = c, c>. 



(30) 



We now have the U-U equation in a form that allows us to 
determine the microscopic frequencies as a perturbation 
expansion in powers of k. However, before we are able 
to determine the transport coefficients, we must relate 
transport coefficients to the frequencies of the linearized 
hydrodynamic equations. 
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where c„ and cp are the specific heats at constant density 
and pressure, respectively, k is the thermal conductivity, 
r\ is the shear viscosity, and £ is the bulk viscosity. The 
speed of sound is given by 
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For the non-condensed dilute boson gas, the specific 
heats are easily obtained from the grand potential 
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From this we derive the the entropy density, 

The specific heats are then given by 
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In the next section, we derive the normal mode fre- 
quencies pip from the linearized U-U equation and 
thereby find microscopic expressions for the transport co- 
efficients. 



V. MICROSCOPIC NORMAL MODE 
FREQUENCIES 



IV. HYDRODYNAMIC NORMAL MODES 

A derivation of the linearized hydrodynamic equations 
and the hydrodynamic normal mode frequencies for a 
fluid of point particles at temperatures above the critical 
temperature for Bose-Einstein condensation can be found 
in [16| (Chapter 8). For such a fluid there are five micro- 
scopically conserved quantities and therefore five hydro- 
dynamic normal modes. The normal mode frequencies 



To obtain microscopic expressions for the hydrody- 
namic normal mode frequencies, we return to the U-U 
equation (|29l) and apply standard Rayleigh-Schroedinger 
perturbation theory, using the wavevector k as the small 
parameter and VTkc z as the perturbation. As a first step 
in applying perturbation theory, we expand lo and \<f>) in 
powers of k so that w = uj^ + few' 1 ' + k 2 io^ + . . . and 
\<p) = |^(°)>+fc|0 (1 ))+fc 2 |0 (2 )) + ... and we then substi- 
tute these equations into (|2T)1) . We then require that the 
coefficients of each power of k vanish separately. 
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A. Zeroth Order U-U Equation 



In the limit when k — 0, (|29l) can be written 



\vp / — — i^wp /• (37) 
From this it is clear that the unperturbed eigenvectors 
are just the eigenvectors of the collision operator 

C. Let us denote the eigenvalues of the collision operator 
as A^ so that 

a(°)\ 



CWf) = \pWp 
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We will denote the five degenerate "zero" eigenvalues 
of C by f3 = 1,...,5. The remaining positive non- 
degenerate eigenvalues of C will be denoted by j3 = 
6, ... ,oo. The five eigenvectors of C corresponding to 
the five "zero" eigenvalues represent the microscopically 
conserved quantities. Properly normalized according to 
Eq. (fl"5j) . they are given by 
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Since these five eigenvectors all have eigenvalues Xp = 
0, we must apply degenerate perturbation theory [l^ll7| 

to determine the particular linear combinations of \4>p°^} 
(j3 = 1,...,5) that will give a well defined perturba- 
tion expansion. Straightforward application of degener- 
ate perturbation theory yields the following appropriate 
zeroth order eigenvectors of C, 
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where a = ^3<7 2 /(5ct ct 4 ). For > 6, |V^° } ) - |^° } ). 

Continuing, we obtain the perturbation expansion for 
the eigenfrequencies of the linearized U-U equation. 

up = - ijXp + VTk(ip^\c„\ip^) 



E 



A/3' - Xp 



(41) 



We now focus on only the frequencies of modes that 
correspond to hydrodynamic modes. To do this, we re- 
strict (3 to only take the values /3 = 1, ... 5. These modes 
represent the conserved quantities and have A^ = 0. The 
restriction in the sum in Eq. (|41l) that Xpi ^ A^ means 
that j3' will only take the values /3' > 6. The sum over 
P' is a sum over only the non-hydrodynamic microscopic 
modes. 

With these considerations, we may rewrite Eq. (1411) as 



up =v T k(i)y\c z \'ij) [ -p > ) 
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From this, one can see that since the eigenvalues \p> 
are all positive, the second order (k 2 ) correction to the 
frequency lies on the negative imaginary axis. Since 
the hydrodynamic modes vary in time as e~ lult — 
e -iw (1) t-|w (2) \t ^ & y\ effects of order k 2 result in damping 
of the hydrodynamic modes towards global equilibrium. 
Thus oj( 2 > represents dissipative effects and will be closely 
related (in fact proportional) to the dissipative transport 
coefficients of viscosity and thermal conductivity. 



B. Calculation of First Order Corrections 

The first order corrections to the hydrodynamic fre- 
quencies are given by 



uf =v T k{i>f\c z \^f). 



(43) 



It is straightforward to compute these quantities in the 
angle basis and we find 
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This result implies the existence of a sound waves trav- 
eling at the speed c s , which matches the result (|3"2"|) ob- 
tained from hydrodynamics. In appendix [S] we discuss 
some of the details of carrying out these calculations in 
the angle basis. 



C. Calculation of Second Order Corrections 

The second order correction to the hydrodynamic fre- 
quencies Up (/3 = 1, ...5) is given by 



,( 2 ) 
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This expression is exact, but requires knowledge of all of 
the eigenvalues and eigenvectors of the collision opera- 
tor. Computation of the eigenvalues and eigenvectors of 
the collision operator was discussed in some detail in Q , 
where it was found that there are an infinite number of 
eigenvalues that converge to an upper limit Am = 03/z. 
This means that the denominator in the sum never be- 
comes large and large groups high order eigenvectors may 
have a cumulative effect on the sum. 

In practice, the eigenvalues and eigenvectors of the col- 
lision operator must be computed numerically. For an 
N x N matrix representation of the collision operator, 
there will be a finite number of eigenvalues that converge 
to the value Am. The remaining eigenvalues will be er- 
roneous because they exceed Am- 

Let us define /3 max as the highest value of j3' for which 
A/3' < Am- We consider the eigenvalues and eigenvectors 
for P' > /3 max to be unreliable due to numerical errors. 

The problem with using Eq. P5|) is that it uses these 
unreliable higher order eigenmodes. This leads to an 
underestimation of the sum, as the eigenvalues A/3' are 
larger than they should be. Still, the contributions of 
these higher order eigenmodes are not negligible, and 
completely neglecting terms with /3' > /3 max will give 
erroneous values. 

We can accomplish this by using a trick that is based 
upon the fact that the eigenvalues converge to the upper 
limit Am- We shall assume that all eigenvalues with j3' > 
Anax are so close to Am that they may be replaced with 
Am in Eq. (|4"1))) . We proceed as follows. 

We can split the sum in Eq. (|4l))) into two parts, 
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where the first term accounts for the eigenmodes which 
are given accurately from the numerical calculation, and 
the second term contains all higher order eigenmodes, 
which are not accurately determined by the numerical 
calculation. 

The second term on the R.H.S of Eq. (|4T)l can be 
further split into three parts, 

00 

E <^i4°Wn^> = 

/S'>/W 

E<4 0) i^i^)^"i^i4 0) } 

/3"=1 
/3»=1 

-E^l^l*)^!^!^)- 
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The first term on the R.H.S. can be identified as 
(■0^|c^|^l ^), due to the fact that the eigenvectors \^a),) 
with ft" = 1, ... 00 form a complete set. 

The second term on the R.H.S can be idendtificd as 



(0)\2 



This follows from the the fact that the 

/,(°)\ 



the eigenvectors \ipa ) were generated from degenerate 
perturbation theory and therefore satisfy {tp^ \ c z \ip^),) = 



SpMi'fl^p fOT /3" = i,..-5. 

Using these simplifications to rewrite Eq. (|4"T|) . we ob- 
tain 
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Am 
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where /j,p> — — . This expression for ujp 

does not depend on any of the unreliable eigenmodes. 

Equation (|49l) is approximate because we have ob- 
tained it by approximating the eigenvalues above /3 max 
by Am- This, however, is a far better approximation that 
using Eq. (l46t , where eigenvalues above /3 max exceed Am ■ 

Even in the extreme case where all eigenvalues and 
eigenvectors are known exactly, we consider Eq. (|4"9")l 
to be preferable because the sum in Eq. (|4"9")l converges 
extremely rapidly. This is due to the fact that the eigen- 
values Api quickly approach Am and thus up/ quickly be- 
comes very large. 

We now divide the calculation of the second order cor- 
rections into two pieces, 

^p = ^f\^f)-^f\c z ^ff, (50) 



and 
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so that 
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Note that this expression does not depend on any nu- 
merically obtained eigenvalues or eigenvectors with /3' > 
/?max, and yet it accurately includes the contributions of 
these terms through Ap. 

Intuitively, A^ represents the result that would be ob- 
tained if all non-zero eigenvalues were assumed equal to 
Am and dp represents corrections due to the fact that the 
first few non-zero eigenvalues are substantially different 
from Am- It is interesting to note that A^ is essentially 
the variance of c z in the hydrodynamic mode \tp^}- 
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It is straightforward to compute values of A^ for 
1, . . . 5 and we obtain 



A 3 
A 4 



2(74 5(To(T4 — 3er| 



= A, = 



12CT2C4 



(53) 
(54) 
(55) 



The calculation of VLp is more tedious but is simplified 
by using the |c, l,m) basis. We now switch to labeling 
the eigenmodes by n,l,m rather than /3'. The range of 
the summation on /?', (6, . . . /3 max ), restricts the range of 
the n summation to values for which < A„./ < Am- We 
then have 



An .1 <Am 
Z.m n,A n />0 



(4 0) i^iCV)«U i^i^ 



1,(0) 



(0) 



>,(0)\ 



(56) 



In carrying out the calculation of for this five cases 
(/3 = 1, ... 5), we find that the only remaining summation 
is over n. The range of n must be restricted so that 
< A nj / < A m- Let us define BmL in analogy with /3 max 
as the largest value of n for which X n j < \ M for a given 
value of We then obtain 



10, 



15(T2 fi n ,2 



I67rg|g ^ x \R n \ 2 

15<T 4 (5cr4O'0 - 3(T§ ) ^ /Lt„,i ' 



(57) 



(58) 



The lower limits of the n summations come from the 
condition that A n ; > 0. More details concerning the 
calculation of A^ and Sip can be found in appendix lAl 
We further define 



Oo = 



\S„ 



47TO-0 

9ct 2 ^-i Mn,0 
n— 2 ' 



(63) 



Note that is identically zero for n > 2. This occurs 

because c 3 ^/ w(c) is a linear combination of t/>oo( c ) and 

■0i°o(c) [see Eq. ([27)]. Thus the function c 3 v /w(c) is 

orthogonal to ^i^C ) f° r n — 2 and the integral in S n 
gives zero for n > 2. Only the two values So and Si are 
non-zero. 

The sum over n in the expression for fio does not in- 
clude the two terms with n = and n = 1. This is 
because Ao,o and A^o are zero. This fact, combined with 
the fact that S n = for n > 2, leads to fio = 0, and as 
we will see, zero bulk viscosity. 



VI. MICROSCOPIC EXPRESSIONS FOR 
TRANSPORT COEFFICIENTS 

By comparing the eigenfrequencies of the U-U equation 
to the hydrodynamic frequencies, we obtain 



n : 



Am 



5fcs(T4(5cr40'o — 3cr|) mnov^ 



Ama 3 , 



A3 
Am 



_ 2mn i4 o 
C, — ii 



(64) 



(65) 



(66) 



$74 — r^5 — 



where 



87TCT0 |i?»| 2 

45 er 4 ^— ' 

n— 1 ' 

ri (2) n (0) 

167rcro \Q n \ 2 4ira a I^nl 2 
45fj 2 ^ AV2 9(t 2 £^ Mn,o ' 



) n = J dcc 3 ^(cj^(c), 





(59) 



(60) 



As discussed in the previous section, O = and there- 
fore the bulk viscosity is identically zero. The fact that 
a monatomic gas has zero bulk viscosity is well known 
[HI, [lH . We have shown that this result follows from the 
conserved quantities and the rotational invariance of the 
collision operator. 

The transport coefficients can be expressed in terms of 
pure functions of the fugacity z by defining 



Then, 



I umkBT 



(67) 



OO 

R n = J dcc*VW) [c 2 - (61) 



S n = / dec 3 v^CoCc) 



(62) 



and 



V _ _&4_ 87TQ-Q y \Q n \ 2 

£ 2cr 3 15z ^ /z„ j2 



"mix 9 

2mK (7(T6<T2 — 5(t|) | 87tcto |-R n | 2 



5fcsC 



4C720"3 



15z fa 



(68) 



(69) 
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FIG. 1: (a) Shear viscosity 77, (b) thermal conductivity re and 
(c) the Euken number versus fugacity in dimensionless units. 
Solid circles indicate the results of Eqs. (|64[) and (|65[) . open 
circles show the first-order classical approximation using one 
Sonine polynomial, and crosses show the results of Eqs. (|64[) 
and ()65[) in the Boltzmann limit z — > 0. 



FIG. 2: (a) Shear viscosity 77, (b) thermal conductivity re and 
(c) the Euken number versus temperature in dimensionless 
units. Solid circles indicate the results of Eqs. (|64[l and (|65fl . 
open circles show the first-order classical approximation using 
one Sonine polynomial, and crosses show the and crosses show 
the results of Eqs. (164|) and (|65[) in the Boltzmann limit z — > 
0. 



VII. NUMERICAL CALCULATION OF 
TRANSPORT COEFFICIENTS 

The numerical calculation of the transport coefficients 
follows closely from a numerical calculation of the eigen- 
values of the collision operator. In Ref. [9|] we discussed 
the method for obtaining the eigenvalues of the U-U 
equation. Of importance is the fact that this calcula- 
tion requires numerical diagonalization of the "collision 
matrix" , which is a finite-sized representation of the col- 
lision operator. All of our results were obtained by per- 
forming the calculation with eight different matrix sizes 
from 490 to 2205, and extrapolating these results to an 
infinite matrix size. 

In Figures [1] and [5] we show our numerical results for 
the transport coefficients. Fig. [T] shows the dependence 



of the dimensionless transport coefficients on the fugacity 
(z) in dimensionless units. Fig. [5] shows the dependence 
on the reduced temperature T/Tq, where Tc is the crit- 
ical temperature for Bose-Einstein condensation. Note 
that the physical transport coefficients have an overall 
dependence on s/T, which is divided out in the dimen- 
sionless transport coefficients. We also plot the dimen- 
sionless ratio mK/{rjc n ), which is known as the "Euken 
number" and is related to the Prandtl number. 



Filled circles depict the "Full" result of Eqs. (|68|) and 
(J5!)|). For comparison, with each curve we have plotted 
two horizontal lines. Crosses depict the "classical Boltz- 
mann" result, which is simply the "full" result at z = 
extended to finite z. Open circles depict the "first or- 
der classical" result that appears in many textbook and 
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is found when the Chapman-Enskog method is applied 
to first order in Sonine polynomials. In the "first order 
classical" calculation, mK/{rjC n ) is exactly equal to 5/2. 

The "first order classical" is noticeably smaller than 
the "full" result. This can be understood by consider- 
ing Eq. (|49|) with M so large (as discussed in Sec. IV Cj) 
that all eigenmodes (even ones that are inaccurate) are 
included. What we find is that as we go to higher order 
eigenmodes, numerical error causes fip' to become neg- 
ative. These negative terms, which are unphysical and 
should not really be included, cause the sum to drop from 
the "full" value to the "first order classical" value. We 
conclude that the "full" value is more accurate and that 
the expansion in Sonine polynomials which produces the 
"first order classical" result underestimates the contribu- 
tions of the higher order eigenmodes. 

The transport coefficients for this degenerate boson gas 
are indistinguishable from the classical Boltzmann values 
temperatures above 20Tc and become sensibly different 
for T < lOTc- As one approaches Tc, the dimensionless 
viscosity decreases to a value that is approximately 82% 
of its value at high temperature. The dimensionless ther- 
mal conductivity shows a minimum around T = 1.5Tc 
at approximately 95% of its high temperature values and 
then increases as the temperature approaches Tc- We 
currently do not have a simply physical explanation for 
this behavior. 



difference to our proper handling of higher order terms, 
which are neglected in the first order classical results. 
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Appendix A: Calculation of A^ and £lp 

In this appendix, we give details of how to obtain Eq. 
(j53")l from Eq. (fST?)) as well as how to obtain Eqs. (|570 . 
([55]) and ([59")) from Eq. (fSTj) . These calculations can be 
done relatively quickly by using the angle basis defined 
in Sec. IIII Bl First we must express the zeroth order 
eigenstates in terms of the |c, I, m) basis. They are 



U7ra w(c) 2 
— ^ c d;,i (<Vi 



/ ; i / 0\ ■ /47ro- w(c) 2 

(c, I, m\ip 2 ) = i J — — c d;,i (d m ,i 



(Al) 



(A2) 



VIII. CONCLUSION 

We have derived an expression for the transport coef- 
ficients of a monatomic dilute Bose gas that obeys the 
Uchling-Uhlenbeck kinetic equation. Our expression re- 
lates the transport coefficients to a spectral decomposi- 
tion of the linearized collision operator. This is in con- 
trast to the standard Chapman-Enskog method which 
uses an expansion in orthogonal polynomials. 

The accuracy of our method relies only of the accuracy 
with which the eigenvalues and eigenvectors of the col- 
lision operator can be calculated. For a gas interacting 
with a contact potential, we can exploit the form of the 
eigenvalue spectrum to obtain a formula that does not 
depend on the unreliable eigenmodes and is influenced 
less by numerical errors. 

We have calculated the transport coefficients of a Bose 
gas interacting with a contact potential for several values 
of the fugacity. Overall, we find the expected decrease of 
the transport coefficients with fugacity, owing to the in- 
creased scattering produced by the Bose enhancement 
factors in the collision integral. However, we observe a 
sharp increase in the thermal conductivity as the temper- 
ature nears Tc . We also find that using the U-U equation 
instead of the Boltzmann equation does not change the 
fact the the bulk viscosity for a monatomic gas is zero. 

Our calculations of the transport coefficients at high 
temperature differ from the first order classical values 
that are found in many textbooks. We attribute this 



{c,l,m\i>D 



(c,l,m\ipl) 



2(72 
5c74S2 



cy/4irw(c) ^c 2 - -j^j 5i,o5m,o, 

(A3) 



' 47T<7oU>(c) 



3(72 



-T^-c 2 Si.o + c5 L1 ) S m , 

5(74 / 

(A4) 



II I / 0\ U-K(T Q w(c) ( [2^ 2 \ 

(c,l,m\tp 5 ) = cW — — I W— c difl + cd it i I d rnfl . 

(A5) 

Expressing the operator c z in the \c,l,m) basis pro- 
vides and elegant and economical way to evaluate A^ 
and We can do this by inserting identity operators 
to obtain 



00 00 
c z — J dc^^ J dc J dc\\c,l,m) 

q i,m q I' ,m' 

x (c, I, m|c z |ci)(ci|c', I', m')(c', I 1 , m'\. 

We then use the definitions (I2U1) and (|3TJ|) to get 

00 



(A6) 



l,m V ,m' 



(A7) 



dnczYr(e,t)Yr {OA) 



(c/'.m'l 
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Performing the angular integration in the bracketed term 
involves Wigner 3j-symbols, but can be simplified to 



,(0 



= / rfccV V \c,l,m) 1 3a 2 x /V*, , 

l=0m=-l y I x x 

X (Ji +1 , m (c,l+l,-m\ + Ji, m (c,l-l,-m\) x / d Cl cly/w(ci)il>n,i(ci) I dc 2C2\Mc 2 )C,z(c2) 







where J hm = yj gpg§+Tj ■ We can also this to get a x _ J mi ^S. 

compact expression for c 2 z in the angle basis. x § l h (S m2 ^t - S„ l2il )5. 

Once we are in possession of the expressions (|A8|) and 
(jAlj) . evaluation of is a trivial exercise. To demon- 



m,m2 

x (^i + l,mi^,Ii+l + Jh,mi^t,h-l) 



strate the method of calculation for Qg, we outline the 2 + 1,7712^2 + 1,1 ,m 2 u l2 — 1,1 ^ ■ 

calculation of Oi below. Starting with Eq. (|5Tj) with (A9) 
/3 = 1, we begin by inserting the expression (IA8|) for 

both occurrences of c z and using the relations in Eqs. Performing the summations and some minor simplihca- 

(|A~T|) and ([26]) to obtain tions, we get Eq. ((57)) . 
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